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Abstract 

We study predictions for the reaction e + e~ — ► 1/1/(717). The complete one-loop 
corrections are taken into account and higher order contributions, in particular those for 
the observed real photons, are added whenever necessary. The event generator /C/CMC, 
a general-purpose Monte Carlo generator for the process e + e~ — > ffwy based on the 
method of exclusive exponentiation, is used as the environment. We extend its appli- 
cability to the process e + e~ — ► uifiwy, I = e,fi,T, where the observation of at least 
a single 7 is required. The exponentiation is implemented in much the same way as 
for the s-channel process alone. In particular, all photonic effects present in the case 
of W exchange, which cannot be included in the s-channel exponentiation scheme, are 
calculated to a finite order only. The real hard photon matrix element is calculated up 
to 0{a 2 ). Leading logarithmic contributions of the two-loop corrections and one-loop 
photonic corrections accompanying real single photon emission are included. The elec- 
troweak corrections are calculated with the DIZET library of the ZFITTER package. 
Numerical tests and predictions for typical observables are presented. 
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1 Introduction 



For the final LEP2 data analysis the total cross-section for the process e + e~ — > vvwy will 
have to be calculated with a precision of 0.5%-l%, and arbitrary differential distributions 
of observable photons also have to be calculated with similar precision In future, for a 
high luminosity linear electron collider like TESLA, the precision requirements will be even 
more demanding. These requirements may be fulfilled by transferring our expertise from 
the case of charged lepton production to that of neutrino production, in particular the more 
involved case of v e v e with t-channel exchanges. The present paper marks an important step 
in this direction. It contains the necessary extensions of the Monte Carlo program /C/CMC 
of ref. [Q], originally written for e + e~ — > //, f — fj,,r, u, d, c, s, b, to the neutrino case. 

In the neutrino pair production process 

e + e~ — > vvwy, (1-1) 

one is interested in observables where at least one high-p r photon is observed; neutrinos 
obviously escape detection. From a methodological point of view, however, it is convenient 
to consider 

e + e~ — > vv (1-2) 

as a (non-observable) Born process, and to incorporate (observable) radiative corrections 
into it, in particular the real photon emissions, which provide the detectable signature. A 
convenient method of exponentiation is discussed in this framework. In order to achieve 
the 0.5% precision level for the vvy final state, the leading-logarithmic (LL) corrections 
have to be calculated up to two or three loops for the virtual corrections and up to two or 
even three hard photons for multiple bremsstrahlung. Mixed real-virtual terms such as the 
loop corrections to real photon emission have to be calculated as wellQ. Needless to say 
that a sufficiently precise integration over the multiphoton phase space within the detector 
acceptance is also necessary. The Monte Carlo (MC) event generator approach is the only 
practicable solution. 

As in the case of any other two-fermion final state in e + e~ scattering, it is possible to de- 
fine certain computational building blocks. Our case does not require the complete two-loop 
effects and we can separate the calculation into two parts: (i) QED: interaction of photons 
with fermions as well as WW'y and WW'yy interactions; and (ii) the rest: non-photonic 
weak and QCD corrections. The type (ii) corrections can be hidden in a few effective cou- 
pling constants. 

The Monte Carlo method is used for the numerical integration over the Lorentz-invariant 
phase space, as usual. The Monte Carlo event generator K.JC MC is documented in ref. O]. 
For a detailed description of its matrix elements for the e + e~ — > ff(rvy) processes we refer 
the reader to refs. 0,^]. 

In section 2 we discuss the implementation of the electroweak corrections. The package 
ZFITTER []5], P is used for this purpose. Basic numerical tests of the code in the absence of 
photonic effects are described. 

'The genuine weak one-loop corrections are sufficient. 
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In section 3 we introduce the photonic matrix elements. We start from the simplest cases 
of and v T production. We then explain the extension of the matrix elements used in the 
CEEX exponentiation of refs. 0-0], to the case of e + e~ — > i> e v e l- The modifications are 
due to the presence of t-channel W exchange. We start with the general description of our 
approximation, and later present the single-photon tree-level amplitude. In particular, we 
explain how single bremsstrahlung amplitudes are used as a building block in the multiple- 
photon amplitudes. Finally we briefly explain the calculation (or construction) of our am- 
plitudes for the different higher order cases: one-virtual and one-real photon, and two-real 
photons. 

In section 4, predictions of /C/C and KORALZ [Q,^]] are given for selected observables of 
the recent LEP MC workshop [p, including results that are used for the final estimate of the 
theoretical and technical errors of our new calculation. 

Section 5 concludes the paper with a statement on the precision of our theoretical predic- 
tions for the uuj process, as compared with the precision targets requested by LEP experi- 
ments [|1|]. 



2 The effective Born approximation for e + e —> vv 

Similarly to the case of pure s-channel two-fermion processes, the electroweak one-loop 
corrections can be incorporated via effective coupling constants of the Z and W to fermions. 
Let us define here the complete electroweak one-loop corrected effective Born cross-section 
for neutrino pair production^. This is by construction a gauge-invariant quantity. It includes 
s-channel Z exchange for all three neutrino species, while for v e pair production it also 
includes t-channel W exchange: 

da- = da ( e+e ~ Zl 3 a + a t + a (2 1) 

The improved Born cross-section originates from a neutral-current matrix element Aiz\\ 

G 

M Z = PeuXz(s) [U^^ (V e + 75) W e ] X (1 + 75) U u ), (2.2) 

and, for u e v e production only, additionally from a charged-current matrix element Mw QH)|]: 

G 

M w = pZ e Xw(t) [M e 7 M (1 + 75) u e ] x [u v j^ (1 + 75) Uu]. (2.3) 



We use here the notations a e = a v = 1, Q e = —1, s w = 1 — M w jM\, and have only three 
form factors K e , p^ v , and p^ v : 



l-A\Q e \s 2 w n e . (2.4) 



2 The notations in this section follow closely those of the package ZFITTER. 
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In the Born approximation, it is p = k = 1. The kinematical invariants are used in the 
approximation m e = 0: 



^1 — cost?), 



u 



We also use: 



Xb{s) 
M%{s) 



-S-t = --(1 + COS0). 



Ml 



-s + M%{s 



Ml-iMn Tn(s) 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



and the width in the s-channel may be chosen constant or s-dependent. The resulting cross- 
section contributions are: 



Or 



O st 



O t 
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\xz{s) p 
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1 + cos ■&) (l + M 2 )+4 cos § 3?e v. 



Me{ Xz (s)x*w(t) Pi PZ*(1 + costf) 2 (1 + v e )}, 



32tt 

sG l L, m^l 2 ^-, n2 

167T 



xw (t) P rj (i+cos^. 



(2.9) 
(2.10) 
(2.11) 



The weak neutral form factors n e and p^ v are discussed in detail in [§, [TT]] and in ref- 
erences therein^. For latest comparisons and applications at LEP1, see also [fTSp, and at 
LEP2, see ref. [JU]. There is one modification with respect to earlier applications, which 
has to be clarified here. Even though the complete virtual corrections are not infrared- finite 
(because of photonic diagrams), we prefer not to split away the photonic part of the virtual 
corrections in our formulae. Instead, we prefer to combine the complete virtual corrections 
with real bremsstrahlung from the package^ fCK, MC. For this reason, we calculate with the 
weak library DIZET of the package ZFITTER the complete virtual correction p^ v . Since 
in the quantity XROK ( 1 ) the QED part is explicitly subtracted, we have to re-establish its 
contribution here according to the formula: 



QEDJJC 



XR0K(1) + QEDJJC 

(A 



In 



(2.12) 
(2.13) 



with L e = \n(s/m 2 e ) and A a finite photon mass. 

The charged current form factor may be extracted from derivations done for ep scatter- 
ing [|T3|-|r5|] . Again, an infrared-finite quantity ROWBf] was constructed, rather than the full 



3 These form factors are calculated in the library DIZET as variables XROK ( 2 ) and XROK ( 1 ) by calling 
subroutine ROKANC (u, - s, t) for s > and t,u < 0. 

4 Of course ICK, MC includes virtual QED corrections of its own. The subtraction of the QED part from the 
complete loop corrections will be defined in section 2. 1 . 

5 With a call on subroutine RHOCC (u, -t, s) (for s > and t,u < 0). This infrared-finite term was 
determined by a subtraction of (gauge-dependent) photonic corrections, which were combined for applications 
at HERA with real-photon corrections. Note also that we use for this part the names of subroutines (but not of 
variables) from the package HECTOR [|0p. 
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form factor p^: 



rev 



QED_CC 



ROWB + QED_CC, 



id; 1 . 2 '<') 



-IL.- -14 + -in 



2 s 



(2.14) 
(2.15) 



For practical reasons, we have decided to insert the neutral current term QED_NC (instead of 
QED_CC) also in the charged current case: 

pZ = ROWB' + QED_NC, (2.16) 
ROWB' = ROWB + 0ED_CC - QED_NC. (2.17) 

All the virtual corrections described here come with ZFITTER v.6.36 (21 June 2001) [@,@]. 



2.1 Form factors and /C/C 



The effects due to the loop diagrams given in formulae Q2.12| ) and Q2.16D are separated into the 
finite parts, which are encapsulated in two of the electroweak form factors and pretabulated 
in /C/C in order to save computer time, and the infrared-divergent part QED_NC which defines 
the (now universal) genuine QED corrections. The term QED_NC is not included in the form 
factors, but enter the QED part of the calculation (encapsulated in /C/C MC). The remaining 
finite parts modify Z and W couplings to fermions as in Refs. 

We are currently not aiming yet at a calculation of the complete second-order corrections. 
We may, therefore, incorporate the difference between the QED parts for s- and t-channel 



JCC-NC 



QED_CC - QED_NC 



—Q 

2tt 



2 s 



+ -ln 2 --4Li 2 (l) + 2 
— s 



(2.18) 



into the electroweak formfactor ROWB', as mentioned above. The difference ( [2.1 8D is infrared 



finite, numerically small, and not enhanced by the large logarithms with respect to f-. The 
numerical contribution is in fact below 0.5% of the Born cross section for LEP2 energies, 
when integrated over neutrino angular variables. 



2.2 Technical tests for the "academic" event selection 

Before the actual discussion of the bremsstrahlung part of the generator, let us make certain 
elementary numerical tests of the implementation of the function of eq. Q2.16| ) within the 
/C/CMC program. This will include the neutral current form factors as well, and will be done 
with the aid of the semi-analytical package /C/Csem, the internal testing program of /C/C MC. 
As a first step we compare the effective Born predictions as calculated by /C/Csem and by 
ZFITTER. In table [I] we present the predictions from the two programs for the processes 
e + e~ — > p + p~, e + e _ — > v^v^, and finally for the processes e + e~ — > v e v e . In all cases the 
agreement between /C/Csem and ZFITTER we find to be sufficiently good, slightly less pre- 
cise in case of e + e~ — ► v e v e . In the latter case possible uncertainties of numerical integration 
(and pretabulation) can be the reason, as the distributions of the scattering angle peakes in 
the forward region. We conclude that the program /C/C MC properly exploits the electroweak 
form factors of the package D I ZET. This opens the way to a more complete treatment of the 
EW corrections in /C/C MC for the neutrino channels. 
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Table 1: Electroweak effects including pretabulation as implemented in /C/C MC. Total cross- 
sections and forward-backward asymmetries are calculated without QED corrections. For every total 
energy the upper entry is from XXsem and the lower one from ZFITTER. 



Channel 


Cms energy [GeV] 


a 


A FB 


e~^e~ — * uf^ ji~ 


91.19 


2.00345818887006 
2.00339325242333 


0.01782302636524 
0.01783385260585 




100.00 


0.05261755867120 
0.05261818913548 


0.58632495479977 
0.58632716129114 




140.00 


0.00697974693325 
0.00697977362831 


0.66477236253951 
0.66477185420337 




189.00 


0.00337662390496 
0.00337666453508 


0.56552465469245 
0.56552686535458 




200.00 


0.00298389425320 
0.00298394853103 


0.55492797790038 
0.55493448086923 




206.00 


0.00279941156110 
0.00279947462003 


0.54984466446977 
0.54985275218873 


E^e~ — ► U,,v,, 


91.19 


3.97432928812849 
3.97416166784039 


0.10951545488039 
0.10954931699922 




100.00 


0.08476335609986 
0.08476366520697 


0.10779125894503 
0.10782386820347 




140.00 


0.00375688589078 
0.00375689763895 


0.10221239752831 
0.10224383875601 




189.00 


0.00115026771057 
0.00115029626002 


0.09199984287991 
0.09203384569328 




200.00 


0.00096201292991 
0.00096202811147 


0.08639737806323 
0.08642343541349 




206.00 


0.00087902320002 
0.00087903346591 


0.08379052580837 
0.08381353397919 


e~^~ e — — * v e v e 


91.19 


3.98484724572530 
3.98468359582307 


0.11158999944974 
0.11162475358080 




100.00 


0.15558776383424 
0.15560448223749 


0.44218368752900 
0.44225943946968 




140.00 


0.04084986326864 
0.04086399163970 


0.82963857483597 
0.82969877450851 




189.00 


0.03937265031133 
0.03939257132943 


0.91661747898646 
0.91665886557827 




200.00 


0.03971937632138 
0.03974060496062 


0.92595698881602 
0.92599592232744 




206.00 


0.03993258937665 
0.03995453232877 


0.93040498709775 
0.93044271454420 



3 Exponentiation and ^-channel W -exchange 

The coherent exclusive exponentiation (CEEX) was introduced in refs. 0,0]. It is deeply 
rooted in the Yennie Frautschi Suura (YFS) exponentiation JT^], and for applications to nar- 
row resonances see also earlier related refs. JT7|, |r"8[]. The exponentiation procedure, i.e. 
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the re-organization of the QED perturbative series such that infrared (IR) divergences are 
summed up to infinite order, is done for both real and virtual emissions at the spin-amplitude 
level (the case of CEEX), while the actual cancellation of the IR divergences always occurs 
at the integrated cross-section level. CEEX is an extension of the traditional YFS expo- 
nentiation, in the sense that, in the standard YFS exponentiation (which we call EEX - for 
exclusive exponentiation), the isolation of the real IR divergences is done after squaring and 
spin-summing spin amplitudes, while in CEEX it is done before. In the actual implementa- 
tion of the CEEX, all spin amplitudes for the fermion pair production in electron-positron 
scattering are handled with help of the powerful Weyl spinor (WS) techniques. There are 
several variants of the WS techniques. We have chosen the method of Kleiss and Stirling 
(KS) [[19|,|20|], because we found KS method well suited for constructing the multiphoton 
spin amplitudes. In refs. 0,0] more details of the approach are available. In particular we 
take all notations and definitions from these works. In the following we recall only the very 
basic formulae of refs. [[3[J|], before we show how the multiphoton CEEX spin amplitudes 
for the W contribution t-channel is constructed and calculated. 



3.1 The master formula 

Defining the Lorentz-invariant phase space as 

J dU V s n {P;p x ,p 2 ,...,p n ) = J (27r)^-E^n (^)%o' 

we write the general CEEX total cross-section for the process 

e + (p a ) + e~(p b ) -> f{p c ) + f{p d ) + 7 (A;i) + 7 (fc 2 ) + ... + j{k n ), n = 0, 1, 2, 00, 

(3.2) 

with polarized beams and decays of unstable final fermions being sensitive to fermion spin 
polarizations (neutrinos are, of course, taken as stable), as follows [Bp: 



where 



fmx( 



~H\ / dLi P S n+2 (Pa +Pb\Pc,Pd, fa,..., fa) Pcl EX {Pa,Pb,Pc,Pd, fa,..., fa) 

(3.3) 



p ( cLx(Pa,Pb,Pc,P d , fa, fa,..., fa) = ±- e Y^ Pa ,.., Pd ) Q (n) J- J2 



3 



ErV'ff' - n j mr( r ) ^ fc l fc 2 k„\ rmj(r) (pk 1 k 2 k ")]* n L n ™ h l h rn 
t a t b °\aXa U X b X b J JV ri \\a 1 a 2 ■ ■ ■ a n ) l JJl n \\a 1 a 2 ' ' ' a n J 1 ° X c \ c ° X d \ d n c n d ■ 

i,j,l,m=0 

(3.4) 

Assuming the dominance of the s-channel exchanges, including resonances, we define the 
complete set of spin amplitudes for the emission of n photons in 0(a r )cEEX if = 0, 1, 2) as 
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follows: 

n 

pe{/,F}' 1 i=l 



m>(2) /pfci fe n A 
-"^n VAo-i ■ ■ • a„) 

n ( n flW (v k i- Y \ (v k 3h. y \ 

e n^r' k'G^+E — :r pJ + e pj 



pe{/,F}" i=l I i=l s [j] l<j<l<n 5 [j] s [i] 



(3.5) 



In the following, we will explain only those aspects of the above formulae that are new for 
the t-channel W -boson implementations (e.g. the phase space and general form of terms will 
remain untouched); for all the rest we refer the reader to refs. [|3|,0]]. We will start from the 
simplest case and consecutively explain how to add more complicated terms. 

Obviously, since only initial-state bremsstrahlung contributes, sums over the partitions 
drops out. The emissions from W bosons (if present) will not be treated as an additional 
source of emission requiring a special partition for itself. It will play the role of a correction 
to the (infrared- finite) /3-functions at any perturbative order. 



3.2 Case of and v T production 

For all neutrino flavours except u e , we can limit ourselves to changes of the numerical values 
of the Z couplings from quarks or leptons to neutrinos. All the formulae remain unchanged 
otherwise, and the previously estimated precisions for the muon channel remain valid also 
for the neutrino case. The only new element is that of the implementation of electroweak 
form factors for the neutrinos. They were not investigated for the neutrino channel until the 
present paper, for DIZET being either a weak library of /C/C MC or a stand-alone code. 



3.3 The v e implementation, general aims 

In the case of calculations performed at any fixed order, there is, in principle, no need to 
worry about gauge invariance, cancellation of infrared or ultraviolet singularities, etc. How- 
ever, in practice it is sometimes quite non-trivial to achieve these essentials. Already the 
introduction of the W and Z boson propagators require summation of infinite series of par- 
tial contributions from any order of the perturbation expansion. Otherwise, the cross-sections 
at the peak of e.g. the Z resonance would not be well defined. There are standard techniques, 
such as the renormalization group, structure functions, and exponentiation, for summing up 
leading higher-order terms arising from ultraviolet, infrared, or collinear singularities. 

Our ambitions here are rather limited. We want to exploit the relative similarity of the 
spin amplitudes involving t-channel W exchange and s-channel Z exchange at relatively low 
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energies. Even though diagrams for electron neutrino pair production involving W exchange 
with photon lines attached only to electron lines are not gauge-invariant (contrasting with the 
case with the analogous s-channel Z exchange), the necessary contribution for completing 
the gauge-invariant amplitude is small (and even vanishes within the LL approximation). W- 
exchange also drops out at sufficiently low energies, where it is legitimately approximated 
as a contact interaction. As a first step, we shall exploit the simplified amplitudes for the 
W^-exchange in the contact interaction approximation, which formally will look exactly the 
same as the contribution of an additional heavy Z' . As a second step we shall calculate the 
difference of the correct/complete perturbative results at fixed, first and second order, with 
the above approximation. In such a two-step procedure we can easily use the already devel- 
oped and tested formulation for the s-channel CEEX and the corresponding program code. 
In addition, we shall also be capable dividing the VT-exchange amplitude into individually 
gauge-invariant parts of a well-defined physical origin. 

In the following we will explain our approach in more detail. We start with a common 
technical trick exploiting Fierz transformation: for massless neutrinos it is possible to rear- 
range the lowest order e + e~ — > vv H^-exchange amplitude into a form identical to the one 
for Z' production exchange, except that the propagator now depends on t and the coupling 
constants are redefined. The complete Born-level spin amplitude then reads 



03 (I; X) = 03 X) = 23 (X) = ® MM (X) = 

= ie 2 ^ U b( X ) ( G f,p)[ba] (Gf,u)[cd] H B = ^ ^fbc][cd}( X ) '> 
B=-y,Z,W B=^,Z,W 

( g *Jm = v(p b , \ b )G^u{p a , A ), (Gf^)[cd] = u(p c , X c )Gf^v(p d , X d ), 

G£. = 7„E w *0? e > G f,, = ^J2^9x' f , ^a = ^(1 + A 75 ), (3.6) 



A=± A=± 



n B u = w ( x ) - , u 2- 

t — M w 



X 2 - M B 2 + iT B X 2 /M B 



Similarly to the case of Z and 7 exchanges, described in ref. M\, the W contribution takes 
the following form: 

r [ W,e WJ rp rpl , W,e WJ jp tt ] 

<rgW , x \ _ Qjc 2 A Aa 9 ~ Xa cXa bXd Xa SXa XcXb AaAd -l (3 7) 

[feo][cd]v I £ _ ' \ ■ ) 

where t is calculated uniquely for every individual event from the 4-momenta oiQ e + , e~ , u, v, 

W,e 1 n A W,v e _ n 1 

9\=-l,l ~ _ V2sin0 w ' U anCl 9x=-l,l - U ' V2sinfl w - 

Let us now consider the amplitudes involving the emission of real photons. We also apply 
the Fierz transformation. The contribution Q3 147 is then added at any place where the standard 
Z contribution (eq. (44) of ref. [Q]) occurs - that is, in the definition of all f3^\ Note that 
such an approximation can be used at any order of the perturbation expansion. This approxi- 
mation preserves gauge invariance, because the resulting spin amplitudes look formally as a 

6 Except some contributions to f3\ fa, where the "correct" transfers involving momenta of the photons has 
to be re-established. See discussion below. 
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contribution from an additional heavy Z' , with the appropriately chosen coupling constants 
and propagators. 

The recipe for extending the /C/C MC amplitudes to the z/ e -pair production at this in- 
troductory level is relatively simple: modify eq. (43) of ref. [Q] and use it later on as a 
building block for all other amplitudes, virtual corrections and hard bremsstrahlung alike[|. 
For the processes at low energies (substantially lower than 80 GeV CMS energy) such an 
approximation coincides with the standard approximation of the contact interaction, where 
the momentum transfer is completely neglected in the W propagators and the W-'j vertices 
are neglected. 

At higher energies, we have to take the t-channel transfer into account; this modifies the 
structure of the amplitudes significantly. We may still, as an intermediate solution, introduce 
an auxiliary single "mean" ("effective") transfer t for the entire amplitude (including prop- 
agators in loop diagrams). With such an auxiliary transfer, the structure of the whole set of 
spin amplitudes still coincides with the one for s-channel processes, provided that we also 
drop all W-j interactions. The introduction of such an auxiliary transfer t is the source of 
certain ambiguity in the case of an event with hard photons. The optimal choice of t should, 
of course, minimize the unaccounted higher-order corrections. 

At present, the routine KinLib_ThetaD of /C/C MC is used to define the transfer t , 
the same as for the calculation of the ^-dependent box corrections in the earlier published 
versions of the /C/C MC. 

3.3.1 One real photon 

The starting point is the well-known 0(a) spin amplitude for the single -photon brems- 
strahlung. We have to consider it anew, because we need it in conventions of ref. [Qj. In 
particular we need to keep track of the relative complex phases of parts of the amplitude, 
which enter the soft photon factors and the remaining finite parts. Also, we want to (re)use 
the part of the amplitude for channel in the v e case without any modifications. This will 
be a starting point for obtaining (3^ of eq. ( |3.5D , which will be later incorporated into our 
general scheme of exponentiation, exactly as explained in [Q] . 

The first-order matrix element from the Feynman diagrams depicted in fig. [1] reads 



7 In the program this is realized by calling subroutine GPS_BornWPlus from subroutine GPS_BornPlus. 



KS =eQ e v(p b , \ b ) M\ d I} Fa \ £ u(p a ,\ a ) 




+eQ e v(p b: X b ) CM Pb ou Wi} u ^ A «) 



(3.8) 
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fWVW 7 



(D 



VVWA/ 1 



CD 




Figure 1: 77ze Feynman diagrams for e + e — > u e u e ^. 



or, equivalently, as follows: 

(£) = -M° + -M w 7 + M 2 WWl + M 3 W7 
M° = eQ e TJ(p 6 , A b ) Mj d /} * a *™~ a h u( Pa , A.) 

+ eQ e i;(p 6 , ^) 4(^0 ~*»t, m+ ft Mft u(p , A.) 

M 1 = +e u(p b , A 6 ) Mm " A a )e* 1 (fci) • (p c -p a +Pb~ Pd) 



t a -M£rt h - M\ 



M* = +ev(p b ,\ b )gYj d v(p d ,\ d )u(p c ,\ c )gZj a h u(p a , A.) 



2 ' 
1 



f a — M-ji tfj — M, 



= -e tJ(p5, A b )^X ft v(p d: X d )u(PcA c )gZZC(ki) «(p«, A a ) 



Wei/ 



1 



1 



2 ' 



(3.9) 



where 



M m = ^ 2 £ n B W^(G^) M 

B=W,Z 

is the annihilation scattering spinor matrix, including final-state spinors and 



(3.10) 



Wev 
j A c ,A a 



l 



sin 9\ 



(3.11) 



7 w 



For the W contribution, the subscripts in M{/} define the momentum transfer in the W 
propagator n^(X): for ac the transfer is t a = (p a — p c ) 2 , for bd it is t b = (pb — Pd) 2 - If both 
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are explicitly marked, then the amplitude 



M^f = ze 2 (G^) M (GW") M (3.12) 



has to be used. 



We split the above expression into soft IR parts proportional to m) and non-IR parts 
proportional to Employing the completeness relations of eq. (A14) from ref. [U to those 
parts we obtain: 

Tyr irn (Pkl) _ _ X^*B \PbPa] TJ \Pak lP al , e Qe \T^y lp h k lPh l jg tp b Pal 

JVL i{/} U<xJ - 2k 1 p a ^ ^pJ [ ] Lpa<xiA a J + ^ LAfcO-ipJ -° UaJ M 

p p 



(3.13) 



The term Aiww-y corresponds to the last three lines[|of eq. ([3.8[). These contributions are 
also IR-finite. At this stage we keep transfers in the t-channel propagators, which depend 
on the way how the photon is attached to the fermion line. The summations in the first two 
terms get eliminated by the diagonality property of U and V (see also ref. [^p) and we obtain 



p 

B ( P kl\ _ _ TPiPal rj [PafclPal , e Q e ST^ y rPtfclPb] m \PbPa] . „ 

r {I} [\ai) - 2k lPa ' " ^AbPaJ [cd\ U [ Pa <Jt\ a \ + 2k lPb ^ U^ipJ 25 IPbXj M 



(3.14) 

The soft part is now clearly separated from the remaining non-IR part, necessary for the 
CEEX. In 23 [ p x ] we use an auxiliary fixed transfer to, independent of the place where the 
photon is attached to the fermion line. In 23 we provide the residual contribution calculated 
as a difference of the expression calculated with the true t-transfers t a , % and the auxiliary 
one t , common to all parts of the amplitude. Note that 23 = 23 + 25. 



8 The term M \ VWl + Mypw-y + M% /Wl originates from the WW-f vertex 

-ie[g^(p - q) P + 9v P {q - r)^ + g vp (r - p) v ] 

where all momenta are outcoming, and indices on outgoing lines are paired with momenta as pP, q v r p ; M\ VWl 
originates from the term where g^ v connects the e~ 



12 



3.4 One- and two-loop and one-loop-one-real-photon QED corrections 



In both the one- and two-loop virtual corrections we use the same formulae as for the 
s-channel. This is a very convenient approximation, because it simply requires the W- 
exchange amplitudes to be multipied by already known functions. As a consequence, the 

(1 2) 

YFS form factor as well as the W contributions to the ' ' functions are readily available 
within the fCKMC environment. 

(1 2) 

The one-loop, complete, virtual electroweak W^-exchange contributions to the (3q ' func- 
tions are known. They are given by the difference between the exact contribution and the one 
for the s-channel, which is given in eq. ( [2.1 8[ ). The above virtual IV-exchange one-loop con- 

(1 2) 

tribution to /3q does not include any numerically sizeable terms with respect to the scale 



of the term defined in eq. fl2.18| ). In the program, the one-loop P^-exchange contribution to 

(1 2) 

(3q ' is located in what we call the electroweak W form factor. 



Encouraged by the smallness of the above IV-exchange one-loop virtual contribution to 

o ' 



(1 2) 

Pq , and inspired by the contact interaction approximation, we assume that the same is true 



(1 2) 

for two-loop- and one-loop-one-real-photon QED corrections in the complete /3q ' . At this 

i 1 



(1 2) 

stage we do not have the complete calculation for (3q . We assume that the approximate 



(1 2) 

[3q ' discussed above and used in this work differ from the complete ones by numerically 
small 0(a 2 ) terms; see also jJZl]]. We will thus use the same s-channel one-loop contribution 

(1 2) 

to Po ' for the W -exchange, amplitudes. This also makes sense because important leading- 
logarithmic photonic corrections are universal, i.e. the same for any hard process. 

The above approximation is, of course, our main source of theoretical uncertainty, which 
we estimate to be of the order of 1% of the cross-sections of the single-photon observables. 
As a guide in estimating the size of the above uncertainty for a given uu{n^) final state, we 
use conservatively the entire size of the one-loop term defined in eq. (|2.18|). 



3.5 Double-photon matrix elements 

Complete double-bremsstrahlung spin amplitudes are at present included. Their contribu- 
tion, as we can see later, turns out to be rather small. We will discuss these exact two-photon 
amplitudes elsewhere, as well as the related questions of the numerical stability; see the 
similar discussion for the single-bremsstrahlung amplitude described in the next section. At 
present, the extensive tests of the type we performed for the s-channel amplitudes, e.g. with 



the calculation based on ref. [|22|], are not yet completed. Gauge invariance was used as 
a main test, so far. We have noticed numerical stability problems, we could not use the 
complete amplitudes in cases when the photon transverse momenta were smaller than some 
fraction of the electron mass; however, we have checked (for a statistics of 800,000 events) 
that this effect is of no numerical relevance. 



3.6 Implementation of photon emission for VK-exchange in KK MC 



Let us now inspect formula Q3.14| ) (see also fig. p. It can be divided into three separately 



gauge-invariant parts. The first gauge invariant part is formed out of the contributions from 
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diagrams (1) and (2) and diagrams (3) and (4) with the common effective transfer^ to: 
s$(ki)*B [ p x ] + r^jy The second gauge-invariant part includes the contributions from di- 
agrams (3) and (4), responsible for the restoration of true t a and t b transfers in place of 
t , combined with part of diagram (5). It also includes the expression rfn + M] VWl of 



eq. d3.14p . The third gauge-invariant part is formed out of the remaining two expressions 



The formula of eq. fl3.14h requires an additional refinement in cases when more than 
one photon are present in the event. In such a case the contribution + AA\ VWl would 
lose gauge invariance because of a "non conservation" of the four momentaF] p a + p b ^ 
p c + Pd + k±. This is cured as follows: if the momentum carried by other photons is in 
the same hemisphere as p a , we choose t a = (p b — k\ — p d ) 2 , t b = (p b — p d ) 2 and e t = 
2e£*(k 1 ) ■ (p b - Pd)», otherwise we assign t a = (p a - p c ) 2 , t b = (p a - k x - p c ) 2 and 
e t = 2e£'*(ki) ■ (p c — Pa)^. In addition, the expression Ai\y Wy is modified as follows: 

UJ " + [XbpJ \t a - M 2 w )(t b - M 2 w y (3 - 15) 

which coincides with the original expression, if additional photons are absent. Also, if the 
additional photons are collinear with beams, then this choice is consistent with including 
them into an "effective beam", in agreement with the principles of the leading-logarithmic 
approximation. 



3.7 Additional pair correction 

The effects from the emission of real charged pairs accompanying the uv'y final states may 
be especially important, because they may change cross sections and distributions through 
the non-trivial experimental event selection (cuts). It should be kept in mind that the typical 
experimental event selection for the neutrino final state requires that there be no charged track 
in the detector (veto). Virtual corrections due to fermion loops in the vertex functions are 
included in the /C/C MC, and can be switched on as explained in the program documentation. 
The corresponding real pair emission then needs to be added also. This can be done by means 
of a separate Monte Carlo generation using any massive four-fermion MC event generator, 
for example KORALW p3|]. 



3.8 Numerical tests using semi-analytical calculations 

Once all necessary ingredients of the /C/C Monte Carlo program are explained, let us start 
numerical studies, first for the inclusive quantities, that is for the integrated cross sections. 

9 In the program, this part of the amplitude is calculated in subroutine GPS_HiniPlus and in 
GPS_HiniPlusW; variables Csuml, Csum2. The contributions yrfn + -M-ww-y) an ^ v^WWi + 
M^), defined later in the text, are calculated as variables Csum3 and Csum4 of subroutine 

GPS_HiniPlusW. 

10 For pure s-channel amplitudes this complication is absent; the propagators of the internal vector bosons do 
not depend on the individual photon momenta, rather on the momenta of the external fermions only, see ref. |Q| 
for details. 
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In table ^ we show comparisons of results from three programs: /C/C MC, /C/C sem and 
ZFITTER. We see good agreement in all three final states, v^ T and //. Only in the case of 
v e ivmax = 0.90; 0.99) was the agreement less satisfactory. These comparisons provide an 
important technical test of our program, although these observables are, of course, academic; 
their definition requires cuts and tagging of invisible neutrinos. 



Table 2: The comparison of three programs: /C/CMC, /C/C sem, and ZFITTER. 



/ 


(a) /C/Csem 


(b) CEEX2/tl (c) ZFITTER 6.36 


(b-a)/a 


(c-a)/a 




<r(w)tpb], «max = 0.01, 189GeV 


V e 

V U 
V r 


25.2091 ±0.0000 
0.7377 ±0.0000 
0.7377 ±0.0000 


25.2351 ±0.0142 
0.7383 ±0.0018 
0.7384 ±0.0018 


25.2127 ±0.0000 
0.7376 ± 0.0000 
0.7376 ± 0.0000 


0.0010 ±0.0006 
0.0008 ± 0.0024 
0.0009 ± 0.0024 


0.0001 
-0.0001 
-0.0001 




2.1645 ±0.0000 


2.1635 ±0.0024 


2.1682 ±0.0000 


-0.0005 ± 0.0011 


0.0017 




ff(«max)tpb], «max = 0.10, 189GeV 


V e 

v a 


32.5172 ±0.0000 
0.9671 ±0.0000 
0.9671 ±0.0000 


32.5312 ±0.0153 
0.9674 ±0.0019 
0.9678 ±0.0019 


32.5572 ± 0.0000 
0.9669 ±0.0000 
0.9669 ±0.0000 


0.0004 ±0.0005 
0.0004 ± 0.0020 
0.0007 ± 0.0020 


0.0012 
-0.0002 
-0.0002 


n 


2.8227 ±0.0000 


2.8234 ± 0.0027 


2.8276 ±0.0000 


0.0003 ± 0.0009 


0.0017 




<T(«max)tpb], «max = 0.80, 189GeV 


v e 
v n 

V T 


45.5520 ±0.0000 
7.6919 ±0.0000 
7.6919 ±0.0000 


45.4818 ±0.0166 
7.7001 ±0.0039 
7.7018 ±0.0039 


45.6423 ± 0.0000 
7.6907 ±0.0000 
7.6907 ±0.0000 


-0.0015 ±0.0004 
0.0011 ± 0.0005 
0.0013 ± 0.0005 


0.0020 
-0.0002 
-0.0002 


[J, 


6.9623 ±0.0000 


6.9658 ±0.0031 


6.9740 ± 0.0000 


0.0005 ± 0.0005 


0.0017 




<r(iw)[pb], « max = 0.90, 189GeV 


V e 
V T 


45.6631 ±0.0000 
7.9012 ±0.0000 
7.9012 ±0.0000 


45.6038 ±0.0166 
7.9100 ±0.0039 
7.9117 ±0.0039 


46.1404 ± 0.0000 
7.9008 ±0.0000 
7.9008 ±0.0000 


-0.0013 ±0.0004 
0.0011 ± 0.0005 
0.0013 ± 0.0005 


0.0105 
-0.0001 
-0.0001 




7.2070 ± 0.0000 


7.2115 ± 0.0032 


7.2192 ± 0.0000 


0.0006 ± 0.0004 


0.0017 




CT(« max )[pb], « max = 0.99, 189GeV 


V e 

V f. 
V T 


45.6770 ± 0.0000 
7.9102 ±0.0000 
7.9102 ±0.0000 


45.6247 ±0.0166 
7.9190 ±0.0039 
7.9207 ±0.0039 


46.3844 ± 0.0000 
7.9100 ±0.0000 
7.9100 ±0.0000 


-0.0011 ±0.0004 
0.0011 ± 0.0005 
0.0013 ± 0.0005 


0.0155 
-0.0000 
-0.0000 


(1 


7.6156 ±0.0000 


7.6234 ±0.0032 


7.6282 ±0.0000 


0.0010 ± 0.0004 


0.0017 



4 Numerical results from /C/C MC and KORALZ 

Let us now present new important numerical results using observables (event selections) 
defined in the LEP2 MC Workshop [[[]], without any modifications. These observables were 
defined to suit the needs of all LEP experiments and to match their particular experimental 
studies. 

Observables can be divided into two groups: single photon tagged observables Nul, 
Null, Nul2, Nul3, Nu3, Nu4g, Nu7, where the precision target required by experiments 
is 0.5%, and the observables with two observed photons Nu2, Nul4, Nu5, Nu6, Nu8, Nu9, 
NulO, where the precision target is 2%, see ref. [jj]]. 

Numerical results are collected in Table 3. For every observable and a total energy of 
189 GeV, numerical results are collected from the /C/C MC with the exact double hard photon 
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matrix element (CEEX2) and the exact single matrix element only (CEEX1). First, we can 
conclude that the predictions from KORALZ presented during the LEP2 MC Workshop and 
for all neutrino observables were generally correct within 4%, in agreement with what has 
been declared. One can even get an impression that the agreement is better, of the order 
of 2%. This, however, is misleading, since these small differences were due to accidental 
cancellations of several effects, all of them of the order of 2-3% . 

Let us now discus results from the /C/C MC with the various options for the matrix ele- 
ment, having in mind an estimate of the theoretical uncertainties. We have found for all ob- 
servables that the correction from additional virtual fermion pairs was systematically —0.6%. 
(The corresponding entries are not included in the actual version of table 3.) Since the virtual 
fermion pairs contribution is rather small, in comparison with other uncertainties, we will use 
the size of the above virtual correction of 0.6% as a conservative estimate of the correspond- 
ing systematic error. Concerning the missing higher order QED contributions (see section 
3.4) we estimate it to be 1% for v e and 0.5% for other neutrinos. Finally we take 0.4% as for 
the systematic uncertainty of electroweak corrections of ZFITTER. The final uncertainty for 
the single-photon observables is thus estimated to be 1.3% for electron neutrino final states 
and about 0.8% for the other neutrino species. 

The uncertainty for observables with two tagged hard photons is higher, probably of the 
order ofQ 5%. As a numerical input we use standard initialization parameters of the /C/C MC 
defined in ref . [|J . 



5 Conclusions 

We have extended the Monte Carlo program /C/C to the neutrino mode. The systematic error 
is estimated to be 1.3% for v e v e ^ and 0.8% for v^v^ and v T v T ^. For observables with 
two observed photons we estimate the uncertainty to be about 5%. These new improved 
results were obtained thanks to the inclusion of non-photonic electroweak corrections of the 
ZFITTER package and due to newly constructed, exact, single and double emission 
photon amplitudes in the /C/C MC for the contribution with the t-channel W exchange. The 
virtual corrections for the W exchange are at present introduced in the approximated form. 
The exponentiation scheme CEEX is the same as in the original /C/C MC program of refs. 0- 

Let us also note that in the other Monte Carlo programs available in the literature for the 
neutrino channel, see eg. refs. [ |2"4"| , |2~5| ], the exact double-photon bremsstrahlung amplitudes 



are also included. Direct comparisons with these programs should be done at a certain point. 
This could provide an independent cross check of our double-photon matrix element and 
may also shed light on certain non-negligible differences between the results for the neutrino 
channel LEP observables collected from various MC event generators during the LEP2 MC 
Workshop [JTp . Therefore, a better answer could be provided on the total theoretical and 
technical uncertainties in these calculations. 



1 1 Such an estimate was obtained by switching on and off the double-photon contribution to the matrix ele- 
ment. See entries CEEX2 and CEEX1 in Table 3. Also, the different options of reduction procedures could 
be compared and the corresponding contribution to the net uncertainties could be reduced. We leave, however, 
these points for future studies. 
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Table 3: Numerical predictions for observables of the LEP2 MC Workshop at a total energy of 
189 GeV (input data of the Workshop). For each observable results are shown in 3 lines: (i) CEEX2 
0(a 2 ) KK MC, (ii) CEEX1 0(a r ) KK MC, (iii) KORALZ MC. The difference 6 is shown in the 
last column: it is the deviation from 1 (x 100) of the ratio of a given result to that of CEEX2. 



Label of obs. 


Program 


Cross section [pb] 


5 ratio 


Nul 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


3.1710-10° ±7.43-10~ 3 
3.2139-10° ±3.85-10~ 3 
3.2244-10° ±4.34-10~ 3 


1.35 
1.68 


Nu2 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


2.1886-10 1 ±9.26-10~ 4 
2.1655-10~ 1 ±8.79-10 -4 
2.1733-10" 1 ±1.20-10" 3 


-1.06 
-0.70 


Null 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


9.1551-10- 1 ±2.26-10" 3 
9.1066-10 -1 ±2.23-10~ 3 
9.1767-10- 1 ±2.42-10" 3 


-0.53 
0.24 


Nul2 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


1.8242-10° ±3.09-10~ 3 
1.8359-10° ±3.07-10~ 3 
1.8442-10° ±3.37-10~ 3 


0.64 
1.09 


Nul3 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


1.7775-10° ±6.95-10~ 3 
1.8085-10 u ±2.82-10 6 
1.8045-10° ±3.35-10~ 3 


1.75 
1.52 


Nul4 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


2.0634-10- 1 ±9.02-10~ 4 
2.0267-10 -1 ±8.40-10~ 4 
2.0121-10" 1 ±1.15-10" 3 


-1.78 
-2.49 


Nu3 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


4.2434-10° ±7.78-10~ 3 
4.2912-10° ±4.48-10~ 3 
4.2885-10° ±4.89-10~ 3 


1.13 
1.06 


Nu5 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


1.2128-10" 1 ±6.77-10~ 4 
1.1912-10 1 ±6.26-10 4 
1.1850-10- 1 ±8.83-10~ 4 


-1.78 
-2.30 


Nu6 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


5.7331-10" 2 ±4.64-10~ 4 
5.5817-10 1 ±4.22-10 4 
5.6543-10" 2 ±6.11-10 4 


-2.64 
-1.37 


Nu7 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


4.4676-10° ±7.82-10~ 3 
4.5206-10° ±4.54-10~ 3 
4.5109-10° ±5.00-10~ 3 


1.19 
0.97 


Nu8 


KKMC CEEX2 
IS. IS. ML CiiiiAl 
KORALZ 4.04 


1.7593-10" 1 ±8.31-10 4 

1. / IDZ'IU ± / .03- 1U 
1.7074-10" 1 ±1.06-10~ 3 


-ZAj 

-2.95 


Nu9 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


7.6434-10- 2 ±5.39-10~ 4 
7.4473-10~ 2 ±4.89-10~ 4 
7.4208-10" 2 ±7.00-10~ 4 


-2.56 
-2.91 


NulO 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


2.5362-10- 1 ±1.00-10- 3 
2.5107-10" 1 ±9.39-10~ 4 
2.5040-10" 1 ±1.28-10" 3 


-1.01 
-1.27 


Nu4g 


KKMC CEEX2 
KKMC CEEX1 
KORALZ 4.04 


1.9091-10° ±6.99-10~ 3 
1.9379-10° ±2.92-10~ 3 
1.9398-10° ±3.46-10~ 3 


1.51 
1.61 
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